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What quantity controls the Coulomb blockade oscillations if the dot-lead conductance is essen- 
tially frequency-dependent ? We argue that it is the ac dissipative conductance at the frequency 
given by the effective charging energy. The latter may be very different from the bare charging 
energy due to the interface-induced capacitance (or inductance). These observations are supported 
by a number of examples, considered from the weak and strong coupling (perturbation theory vs. 
instanton calculus) perspectives. 
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The quantum fluctuations of charge in a Coulomb 
blockade (CB) quantum dot flu 0] have been a subject 
of intense studies (see Refs. ||, |J for a recent review). 
It is widely accepted that the dimensionless (in units of 
Gq = e 2 /{2n%)) conductance, G, of the interface be- 
tween the dot and a lead controls the strength of the 
quantum fluctuations. If G < 1 (weak coupling) the fluc- 
tuations are small (possibly apart from the degeneracy 
points) and usually may be accounted for by the pertur- 
bation theory in G [|[ ||, |6). On the other hand, for G > 1 
(strong coupling) the CB is basically destroyed by the 
fluctuations. The remained weak CB oscillations may be 
described in the instanton approach [Q, |[ [| [HJ [HI (fj||. 

What happens if the conductance of the interface is 
essentially frequency-dependent, G = G(uj) ? Is it the 
dc conductance, G(0), or rather an ac one, G(fl c ), (or 
may be none of them) that governs the strength of the 
quantum fluctuations ? Is there a single parameter that 
divides the regions of applicability of the perturbation 
theory and the instanton calculations ? What energy 
scale controls the thermal fluctuations ? These are the 
questions we address in the present paper fl3fl . 

We argue that it is the dissipative ac conductance, 
5fteG(f2 c ), that controls the quantum fluctuations and 
seamlessly divides between the two methods of descrip- 
tion. The relevant frequency f2 c has the meaning of an 
effective charging energy of the dot, Q c ~ e 2 /(C + C), 
where C is the bare capacitance of the dot and C(uj) 
is the parallel capacitance induced by the dispersive in- 
terface. Frequency il c is the solution of the following 
self-consistent equation: 

G^ (C + <5(fi c )) = (I- 1 ; C(u) = G Q d u G{iw) . (1) 

The effective charging energy, f2 c , serves also as the char- 
acteristic scale for the thermal smearing of the CB. The 
induced capacitance, C(w), is not necessarily positive and 
therefore O c may be on the either side of the bare charg- 
ing energy, E c = e 2 /(2C). 

To illustrate these points we consider three specific ex- 
amples. In one of them (tunnelling barrier): G(0) ~ 
G(tt c ). In the second (resonant impurity interface) the 



dc conductance is parametrically larger than the ac one, 
G(0) ^ G(ft c ). Finally, in the third example (shallow 
2d lead) the situation may be reversed: G(0) <C G(il c ). 
We treat all the examples from both weak and strong 
coupling perspectives. While the former is relatively 
straightforward, the latter requires certain technical ad- 
justments to the way the instanton calculus is usually 
applied to the CB problem. 

The CB is described by the following partition function 

OC r. 

Z(q) = ^ W / V< P > ( 2 ) 

W =~°° 0(/3)-0(O)=2ttW 

where q is the external charge controlled by the gate volt- 
age and W are integer winding numbers The imagi- 
nary time phase action has the form [jl4| 

S[<f>] = Jdr - JJdrdT' e*WK(j - T')e~ 1 ^ . 

o ' o 

(3) 

The boson Matsubara transform K (ui m ) — K m of the 
kernel K{t) is given by 

^ m = Tri^W n g«Wt{ a W m -^)} , (4 ) 

where W n = W(e„) is the tunnelling matrix connect- 
ing the dot and the lead and = Q^ d ^{e n ) are 
the fermion Green functions of the dot and lead corre- 
spondingly. The linear ac conductance of the interface 
is determined by the kernel, K m . Indeed, employing the 
Kubo formula, one finds for the Matsubara transform of 
the conductance G m = —Aiv\oj m \~~ l $leK m . After ana- 
lytical continuation from the upper complex half-plane, 
uj m — > — ioj, one obtains the complex ac interface con- 
ductance G(ui). The kernel K m is also related to the 
non-linear current-voltage characteristic of the interface 
|TJ||. Namely, I(V) = 2e ^m{K \ Um ^ ie v } is the cur- 
rent through the interface, provided that the voltage V 
is applied across it. We shall mostly employ the ac inter- 
pretation, although the I — V language is also possible. 
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In the weak coupling regime one expands exp{— S [(/)}} 
in the power series over the K-term in Eq. (^|) . Perform- 
ing then the Gaussian ^-integration and VF-summation, 
one finds for the average number of electrons on the dot, 
SN{q) = q+ {2/3E c )- 1 d In Z/dq : 



SN(q) 



■DC 

J <IttK{t) 



-E c {\T\-2qT) 



(5) 



We have restricted ourselves to the temperature range 
[3Ec > 1 and thus extended the integration to infinity. 
For a closed dot the electron number stays constant for 
\q\ < 1/2 and jumps by one at the degeneracy points, 
q = ±1/2. A weak coupling, K(t), makes the electron 
number to change smoothly in between the degeneracy 
points. It is clear from Eq. (^|) that, apart from the im- 
mediate vicinity of q = ±1/2 it is the short time, 
|t| < E^ 1 , behavior of K(t) that determines SN(q). 
One expects, thus, Eq. (^) to be valid if fl c f» E c and 
UeG(Sl c ) < 1. 

If this is not the case one may employ the instanton 
approach. To this end one has to find stationary config- 
urations of the action ([|) with a fixed winding number, 
e.g. W = 1. We are unable to execute this program 
for an arbitrary K(t) and employ instead the variational 
method. As a trial variational solution we take a family 
of Korshunov instantons fll7[| : 



i<t>(r) _ 



(6) 



where z is a complex number with \z\ < 1 that determines 
location and the width of the instanton. The action of a 
trial configuration is, in general, complex. Its imaginary 
part is simply a shift of the background charge, q, which 
affects only a phase of the CB oscillations. Hereafter 
we focus on the real part. Substituting the trial phase 
configuration, Eq. (p), into Eq. (0), one finds 



» e S(S!) = |2- 



where /3Q/(2ir) = (1 — | ^; | 2 ) 1 > 1 is the inverse 
width of the instanton. One should minimize now 
Z\ ~ ((3D,) exp{— SRe 5(0)} with respect to the instan- 
ton width, f2 _1 . The factor /3Q originates from the 
integration over the zero mode associated with the in- 
stanton location, argz. The minimum is reached at 
Q = fl c that is a solution of the stationary point equation 
d 3?e S(H)/dSl = At low temperature this equation 

essentially coincides with Eq. (|l|). The optimal instan- 
tons are usually narrow, /3J7 C 3> 1, that justifies the non- 
interacting instanton gas approximation. Summing up 
the gas of optimal instantons and anti-instantons |l2j , 
one finds for the average number of electrons 



6N(q) 



C E~ e 



-^s(n c ) sin ( 27rg + <^) ; ( 8 ) 



with a model dependent proportionality coefficient, c, 
and phase, tp. This result is applicable if 3?e5(fi c ) ^> 1. 

We have also employed the trial configuration of the 
form 4>(t) = 4 arctan (expjfir}), that is known to be a 
solution for the Josephson tunnelling problem |lS[ ]. The 
results differ from those obtained with Eq. (^3j) only by 
numerical coefficients. That makes us believe that our 
variational approach is parametrically accurate. We turn 
now to the applications of these results for the specific 
examples. 

(i) Frequency-independent conductance: G m — Gq. 
This is the case e.g. if both lead and dot have continuous 



spectrum, Gr?'^ 



sign(e„), and energy indepen- 



dent tunnelling matrix elements, W. In the low temper- 
ature limit the kernel has the form K(t) — Gq/ (2ttt) 2 . 
Employing Eq. (0) , one finds for a weakly coupled dot || 
SN(q) = ,^R 2 In • In the strong coupling regime the 
trial phase, Eq. (0), is the exact solution of the saddle 
point equation. For fl <C E c the action $teS(Q) — Gq/2 
is fi-independent. As a result, there is an additional 
zero mode associated with the instanton width. The 
presence of such zero-mode complicates calculations of 
the pre-exponential factor. We quote here only the re- 
sult {| H ||: cfi c in Eq. (§) must be substituted by 
(2tt)~ 1 E c Gv \n((3E c ). The additional zero-mode leads to 
the instanton-instanton interactions that break the non- 
interacting instanton gas approximation at temperature 
ks E c Gq exp{— Go/2}. At smaller temperature the 
amplitudes of the higher g-harmonics are probably of the 
same order as the first one. That may lead to a nonana- 
lytic behavior at the degeneracy points, q — ±1/2, in the 
zero temperature limit, similar to that of Ref. p9[ . As 
expected, the unit conductance, Go ~ 1, separates the 
strong and weak coupling regimes. 

(ii) Resonant impurity. Consider an interface where 
the particle exchange between the dot and bulk takes 
place only through resonant impurity levels. For a single 
channel this model was solved in Ref. (2TJ. We shall as- 
sume for simplicity that all the levels are at the Fermi en- 
ergy and have the same coupling to the bulk (width) , T <C 
E c . The Green function of the lead is given by that of 
the resonant impurity level: Qn = (ie n + *rsign(e„)) _1 . 
With the help of Eq. (0), one finds for the kernel: 



K m = - — — In 

4-7T 



1 



r 



(9) 



Accordingly the linear conductance takes the form 
G(oj) = GoTiuj^ 1 ln(l — iuj/T), where Go is the dc reso- 
nant conductance proportional to the number of impuri- 
ties and their coupling to the dot [|l). The short time, 
Tt < 1, limit of the kernel is K(t) = G r/(87r|r|). Em- 
ploying Eq. (||), one obtains for SN(q) not too close to 
the degeneracy points, q = ±1/2 pfJ: 



SN(q) = 



1 



GqT 

8nE c \l-2q l + 2q 



(10) 
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Notice, that the criterion of the weak coupling is 
GqT/E c -C 1 pOj , while the dc conductance, Go, may 
be large. In such a case Eq. (|l|) is solved by £l c s» E c . 
Since ?R.eG(uj) — ttGqT / (2u>) for lu 3> T, one finds that 
SN(q) ~ 5fteG(fi c ). As a result, the weak coupling ap- 
proximation is valid as long as 5ReG(O c ) -c 1. 

In the strong coupling regime, employing Eq. (Q), one 
finds for the variational action 

G r r , n „ 



± denote the sign of 6q. In the small temperature limit, 
[3eo ;$> 1 one finds for the dissipative conductance 



3ie5(fi) = |J 



2il 



e r/n E 1 (r/fi) 



where E\ is the exponential integral. This action pos- 
sess two stationary points as a function of the in- 
stanton width, f2 . One is the stable minimum at 
tt c = ^G q TE c /(2tt) \n 1/2 (G a E c /r), with the action 
5fte S(fl c ) = iril c /E c . The other is the unstable plateau 
at rri -1 > 1 with the action given by half of the dc con- 
ductance: 5fte 5(0) = Go/2, in agreement with (i). One 
may check that fl c is indeed the solution of Eq. (|l|). Up to 
logarithmical factor [Q, Ke5(f2 c ) ~ 5fte G(f2 c ) and the 
condition 5Re G(Q C ) « \JGqY / E c 1 justifies validity of 
the instanton approach. 

At not very small temperature, the stable minima at 
fl c (narrow instanton) is the only relevant one and the 
electron number is given by Eq. (|^). The contribution of 
the wide instantons with O < |T| is exponentially smaller 
~ G(0) 2 ln(/3|r|) e~ G (°)/ 2 and may be disregarded as long 
as G(0) > UeG(Cl c ) « SJe5(Q c ). Nevertheless, due 
to the phenomena mentioned in (i), the wide instantons 
may prove to be important in an immediate vicinity of 
the degeneracy points and at very small temperature, 
I3- 1 < £; c G^exp{-G(0)/2}. If so, then G(0), rather 
than 5Re G(Q C ), controls the CB in that narrow range. 

On the high temperature side the CB oscillations are 
washed out at T > ft c . Notice, that in the strong cou- 
pling regime fl c > E c and therefore the (weak) CB, 
Eq. (||) , is more "temperature-tolerant" than in the weak 
coupling. The reason is that the induced capacitance, G, 
is negative (it is actually an inductance) and thus the 
effective charging energy is increased. 

(iii) Shallow 2d lead. Consider a quantum dot coupled 
through a tunnelling barrier to a shallow 2d gas. The 
bottom of the 2d conductance band is at energy eo below 
(eq < 0) or above (eo > 0) of the dot's Fermi energy. In 
the former case the dc conductance between the dot and 
2d gas is given by the product of their densities of states 
and the tunnelling matrix elements and is denoted by 
Go- In the latter case the dc conductance is zero. We 
shall examine the effect of the proximity to such a lead 
(eo <C E c hereafter) on the CB in the dot. 

The straightforward calculation yields for the coupling 
kernel: 

K(,\ Ga f 1 ~ AH") e- 6 ^ , e <0; 



(2ttt) 2 L d{T)e~ £aT , e >0. 

The corresponding conductance, G(lo), obeys the symme- 
try relation G_(oj) = Gq — G+(o>), where the subscripts 



5fteG + H 



Go 
~2~ 



M-eo) 1- 



£0 



(13) 



Notice, that KeG(w) = for \u\ < £ if the 2d band 
is empty, eo > 0. On the other hand, in the high fre- 
quency limit, uj ^> |eq 1 , 3?eG(cj) = Go/2 irrespective 
to the sign of eo- For the reactive conductance one ob- 
tains ^smG+(uj) = Go/(27r)/(|cj|/eo), where we denoted 
f(x) =1n\x- 1\ - ln\x + x- 1 \tl\x 2 - 1|. 

In the weak coupling limit, employing Eqs. (j^) and 
(|l2|), one finds for the average number of electrons: 



6N(q) 



_^L/ In^|-ln(l + |e |/ J B c + 2g), e < ; 
(2tt) 2 1 - ln(l + e /E c -2q), e > 

(14) 

(we have subtracted an irrelevant constant 
-(G /(47r) 2 )lnA/ J B c , where A is the 2d band-width). 
For the shallow filled 2d gas (sq < 0) the result is only 
slightly different from (i). Surprisingly, for the empty 2d 
band (eo > 0) one finds a g-dependence of 5N, despite 
of the fact that the dc conductance is strictly zero. This 
dependence is parametrically similar to that of the filled 
2d gas, however does not exhibit singularities, indeed 
\q\ < 1/2. It is associated with the virtual transitions 
between the dot and the empty 2d band. Since eo <C E Cl 
such transitions are possible and the coupling strength 
is determined by Go only. For Go < 1 Eq. ([!]) is solved 
by Q c w E c . The applicability of the weak coupling 
approximation is therefore controlled by the condition 
5JeG(fi c ) = G /2 < 1 (and not G(0) <C 1 !). 

We shall approach now the problem from the strong 
coupling perspective. To this end we employ the vari- 
ational expression, Eq. (R). For eo < the variational 
action is 5Re 5 = G (l - Q/(%e ))/2 for fl < |e |, while 
it approaches die S = G /4 + 7rf2/i? c for f2 > |e |. There 
is a broad minimum at |eo| <C ^ c < E c . As a result, 
!Re5(T2 c ) = Go/4 and the electron number is given by 
Eq. (|8|) (possibly apart from the immediate vicinity of 
q = ±1/2). This result is valid for 5Re G(fi c ) = G /2 > 1. 

For the empty 2d band, Co > 0, the variational action 
is 5ie5 = Gor2/(27re ) for ft <C Co and riches Go/4 at 
Eq ^ Q < E c . Naively, it has a minimum at Q c — > with 
3?e 5(0) = 0. One has to include, however, the fluctuation 
part into the stationary point equation, as discussed after 
Eq. (||). As a result, for Go ^> 1 one finds the stable 
extremum at Q c = 2tteq/Go. The corresponding electron 
number is 5N(q) = q — ceq/ (GqE c ) sin(2Trq). Despite 
being obtained in the instanton approach, the correction 
is not exponentially small. The formal reason is that 
dieS(fl c ) = 1, that makes the above result an order of 
magnitude estimate, at best. Fortunately, one can solve 
the problem in a different way. 

As was shown above the "straight" instantons, 
< / > vk(t) = 2irWr/f3, are the most stable. One can, thus, 
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calculate the partition function as a sum over such instan- 
tons. The action is the Matsubara transform of the ker- 
nel: %leS[<t) W ] = -P3teK{u; w ) « W 2 G /(2(3e ). Em- 
ploying Eq. @, one finds lnZ(q) = -/3g 2 (2^ 2 £ /G ); 
valid at small temperature, /3 _1 <C Eq/Gq. As a result, 



6N(q) = q 1 - 



27T% 
GnE c 



\q\ < 1/2 



(15) 



and there are small jumps by tteq/ (GoE c ) <C 1 at the 
degeneracy points q = ±1/2. The presence of such jumps 
could be anticipated from the absence of singularities in 
the weak coupling expression, Eq. ([l4|). 

Equation ([l5]) states that there is a perfect CB with 
the reduced (screened) charging energy nQ c = 2ir 2 £ /Go. 
Indeed, according to Eq. (Q) the induced capacitance is 
given by C(u>) = e 2 G e /(2iruj) 2 ln(l + (lu/e ) 2 ). As a 
result, one finds C = e 2 Go/(47r 2 £o) > C at small fre- 
quency, lo <C £q. The bare capacitance, C, appears to 
be completely screened and substituted by G. Solving 
Eq. (|IJ) with such G, one obtains f2 c = 2tteq/Gq that jus- 
tifies the small frequency assumption, provided Go S> 1. 
The amazing twist is that KeG(f2 c ) = and quantum 
fluctuations are practically absent. Therefore the larger 
Go is - the weaker is the coupling to the lead ! By increas- 
ing Go one can restore the perfect CB with a very small 
charging energy, however. The (exponential) tempera- 
ture smearing sets in at the temperature ~ Q c E c , 
limiting an indefinite increase of Go- 

Strictly speaking, the system is always in the weak 
coupling regime, SfteG(f2 c ) <C 1 (possibly apart from 
G « 1). The perturbative result, Eq. ( |l~4|), is valid, 
though, only for G « 1. For G > 1 Eq. (|L5|) is to be 
employed. (It is in parametric agreement with the varia- 
tional estimate, although the q-dependence of the latter 
is actually wrong.) The two results match at Go ~ 1. 
Indeed, consider e.g. \ — d8N{q) / dq\ q= o- For Go < 1 
one finds X = (2n 2 y 1 G /{l + e /E c ), while for G > 1 
X = 1 — 27r 2 Eq/(GqE c ). Since £o <C E c , these two (weak 
coupling) expressions perfectly match at Go = 27r 2 . 

To summarize: we have argued that not too close to 
the degeneracy points the quantum fluctuations of charge 
are governed by the dissipative ac conductance 3?e G(f2 c ). 
The effective charging energy, tt c , has to be found as the 
solution of Eq. ([!]). It may substantially deviate from 
the bare charging energy, E c , in either direction. For 
5ieG(f2 c ) <C 1 and low enough temperature, /3tt c 3> 1, 
the charge on the dot is approximately quantized. For- 
mally it stems from the fact that the large winding num- 
bers, W ~ /3fi c 3> 1, give a dominant contribution to the 
partition function. In the opposite limit, 5fte G(f2 c ) 3> 1, 
only the instantons with the smallest VF's are important, 
while Q c acquires the meaning of the optimum inverse 
instanton width. In both cases for P^ 1 > J7 C the thermal 
fluctuations smear the CB. On the other hand, very close 



to the degeneracy points and at small temperature, the 
relevant parameter may prove to be the dc conductance, 
G(0). A detailed study of the last issue is beyond the 
scope of the present letter. 
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